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Abstract. There is robust observational evidence supporting the existence of 5 — 20 Mq 
compact bodies in X-ray binary systems and of 10^ — 10^ Mq bodies at the center of many 
galaxies. All these objects are commonly interpreted as black holes, even is there is no 
direct evidence that they have an event horizon. A fundamental limit for a black hole in 
4-dimensional general relativity is the Kerr bound |a*| < 1, where a* is the spin parameter. 
This is just the condition for the existence of the event horizon. The accretion process can 
spin a black hole up to a* ^ 0.998 and some super-massive objects in galactic nuclei could 
be rapidly rotating black holes with spin parameter close to this limit. However, if these 
super-massive objects are not black holes, the Kerr bound does not hold and the accretion 
process can spin them up to a* > 1. In this paper, I consider compact bodies with non-Kerr 
quadrupole moment. I study the evolution of the spin parameter due to accretion and I find 
its equilibrium value. Future experiments like the gravitational wave detector LISA will be 
able to test if the super-massive objects at the center of galaxies are the black holes predicted 
by general relativity. If they are not black holes, some of them may be super-spinning objects 
with a* > 1. 
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1 Introduction 

Today we have robust observational evidence for the existence of 5 — 20 Mq compact objects 
in X-ray binary systems [1] and of 10^ — 10^ Mq objects at the center of many galaxies [2]. 
The stellar-mass objects in X-ray binary systems are surely too heavy to be neutron or quark 
stars for any reasonable equation of state [3, 4]. At least some of the super-massive objects 
in galactic nuclei are too heavy and compact to be clusters of non-luminous bodies, since the 
cluster lifetime due to evaporation and physical collisions would be shorter than the age of 
the system [5]. All these objects are commonly interpreted as black holes (BHs), since they 
cannot be explained otherwise without introducing new physics. However, there is no direct 
observation evidence that they have an event horizon [6], while there are theoretical argu- 
ments suggesting that the final product of the gravitational collapse may be quite different 
from a classical BH [7-12]. 

In 4-dimensional general relativity, BHs are described by the Kerr solution and are 
completely specified by two parameters: the mass, M, and the spin angular momentum, J. 
Instead of J, one can use the specific spin angular momentum a — J /M ^ or the dimensionless 
spin parameter a* = J/M^ . The fact that these objects have only two degrees of freedom is 
known as "no-hair" theorem [13-15] and implies that all the mass moments, and all the 
current moments, 5/, of the space-time can be written in term of M and J by the following 
simple formula [16, 17]: 

M + i5i = M(^i^y . (1.1) 

As it was put forward by Ryan in [18, 19], by measuring the mass, the spin, and at least one 
more non-trivial moment of the gravitational field of a BH candidate, one over-constrains 
the theory and can test the Kerr BH hypothesis. 
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The possibility of testing the Kerr metric around astrophysical BH candidates with 
future experiments is quite extensively discussed in the recent literature. A constraint on the 
nature of these objects has been recently obtained by considering the mean radiative efficiency 
of AGN [20] . The detection of gravitational waves from the inspiral of a stellar-mass compact 
body into a super-massive object, the so-called extreme mass ratio inspiral (EMRI), can be 
used to put very interesting constraints [21-27]. Since the future gravitational wave detector 
LISA will be able to observe about 10"^ — 10^ gravitational wave cycles emitted by an EMRI 
while the stellar-mass body is in the strong field region of the super-massive object, the 
quadrupole moment of the latter will be measured with a precision at the level of 10~^ — 
10~^ [23]. The metric around super-massive BH candidates can also be probed by observing 
their "shadow" [28-30]. Additional proposals to test the Kerr BH hypothesis involve the 
possible discovery of a stellar-mass BH candidate with a radio pulsar as companion [31], the 
study of the Ka iron lines [32], and the analysis of the X-ray spectrum of a geometrically 
thin and optically thick accretion disk [33]. 

A fundamental limit for a BH in general relativity is the Kerr bound |a>,<| < 1. This is 
just the condition for the existence of the event horizon. The accretion process can spin a 
BH up to a* ^ 0.998 [34] and many super-massive objects in galactic nuclei may be rapidly 
rotating BHs with spin parameter close to this limit [35-37]. Nevertheless, if these objects 
are not the BH predicted by general relativity, the Kerr bound does not hold. Interestingly, 
in this case the accretion process can easily spin them up to a* > 1 [38]. Such a possibility 
is currently ignored in the literature. It is also neglected in those works in which it is not 
assumed that these objects are BHs and it is studied how future observations can test the 
Kerr metric. For experiments like LISA, that rely on matched ffitering, this may be a serious 
problem. 

In Ref. [39-42], I studied some features of the accretion process onto objects with |a*| > 
1. However, an important question to address is if objects with |a*| > 1 can really form. 
In [38], I showed that deviations from the Kerr metric can have the accreting gas spin the 
body up to a* > 1. I used the Manko-Novikov (MN) metric [43], which is a stationary, 
axisymmetric, and asymptotically flat exact solution of the vacuum Einstein's equation. It 
describes the exterior gravitational field of a rotating massive body with arbitrary mass 
multipole moments and has an infinite number of free parameters. However, it has the 
drawback that it is valid only for sub-extreme objects, with spin parameter |a>,<| < 1. So, I 
showed that the equilibrium spin parameter al^ must be larger than 1, but it was impossible 
to discuss the accretion process for a* > 1, compute a*^, and figure out the properties of 
the space-time when a* > 1. An extension to include super-spinning objects, if it exists, is 
non-trivial. In the present paper, I overcome this problem by considering another metric, 
the Manko-Mielke-Sanabria-Gomez (MMS) solution [44, 45]. It is not as general as the MN 
metric, but it can be easily extended to discuss objects with a* > 1. 

The content of this paper is as follows. In Sec. 2, I present the MMS solution adapted to 
the case of super-spinning objects. In Sec. 3, I study the main properties of this space-time 
and I show that the evolution of the spin parameter for an object whose exterior gravitational 
field is described by the MMS metric can be calculated as in the case of a BH. In Sec. 4, I 
discuss the evolution of the spin parameter of an object with non-Kerr quadrupole moment, 
as a consequence of the accretion process. Sec. 5 is devoted to the discussion of the results 
presented in Sec. 4, while summary and conclusions are in Sec. 6. Throughout the paper I 
use units in which Gat = c = 1. 
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2 Exterior field of a rotating body with non-Kerr quadrupole moment 



The Manko-Mielke-Sanabria-Gomez (MMS) metric [44, 45] is a stationary, axisymmetric, 
reflection-symmetric, and asymptotically flat exact solution of the Einstein-Maxwell equa- 
tions. It includes, as special cases, the Kerr and the 5 — 2 Tomimatsu-Sato solutions. It 
is specified by five real parameters: the mass, M, the specific spin angular momentum, 
a — J/M^ the electric charge, Q, and two other parameters, h and ji. The latter of which 
determine the mass quadrupole moment, Q, and the magnetic dipole moment, M. Since 
the MMS solution is reflection-symmetric, all the odd mass moments and the even current 
moments are identically zero. Here I am interested in the vacuum solution only, and I put 
Q — M — which implies /x = 0. The mass quadrupole moment of the gravitational field is 

Q^QKerr-{d-5-ah)M , (2.1) 

where Qxerr — —o?M is the quadrupole moment of a BH and 



M'^-{a- by ' 

In Refs. [44, 45], the metric is written in prolate spheroidal coordinates, which are suit- 
able for slow-rotating objects. It can be adapted to fast-rotating objects by proceeding as for 
the 5 — 2 Tomimatsu-Sato metric [46] : one has to change the prolate spheroidal coordinates 
into oblate spheroidal coordinates. That can be achieved through the transformation: 

X ^ ix , k ^ —ik , (2-3) 

where i is the imaginary unit, i.e. = —1. The line element becomes 



fc2 



+ -{x^ + l)[l-y'')dct>\ (2.4) 



where k — \/—d — 5 and 

The functions A, B, and C can be written in the following compact way [45] 

A = R^ + XiX2S^, B = A + RP + X2ST, C = RT -X^SP . (2.6) 

Here Ai = A;^(x^ -|- 1), A2 = — 1, and 

P = 2kMx[{2kx + Mf - 2y^{2S + ab - h^) - + b^] - 4y^{45d - M^b^) , 
R = 4[fe2(x2 + 1) + S{1 - y2)]2 + (a _ b)[(a - b){d -5)- M^l - y^f , 
S = -4{(a - 6)[fc2(x2 -I- y2) + 2(5y2] + y'^M% , 
T = 8Mb{kx + M)[k^{x^ + 1) + 5{l - y^)] + 

+(1 - y2){(a _ b){M^b^ - 4Sd) - 2M{2kx + M)[(a - b){d -5)- M^]} . (2.7) 
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a* = 0.4 



a* = 1 .4 
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Figure 1. q = {d — S — ah)/o? as a function of the parameter h for a* = 0.4 (left panel) and a* = 1.4 
(right panel). The Kerr solution is recovered for h = ±Va^ — M^. b is given in units of M = 1. 

In general, the mass quadrupole moment of an object depends on its mass and on 
its spin angular momentum in a non-trivial way, according to the specific properties of the 
matter the body is made of. In this paper, I will use the following definition of anomalous 
quadrupole moment q 

Q^-{l^q)a^M . (2.8) 

For (7 = 0, one finds Q — Qxerr, while for g > (g < 0) the object is more oblate (prolate) 
than a BH. Eq. (2.8) is what we expect for neutron stars, with g ^ 1 — 10, depending on 
the equation of state and the mass of the body [47]. It is reasonable to assume that g > — 1, 
since otherwise the effect of the rotation would make the body more and more prolate. For 
given M and a, the MMS solution does not allow for any arbitrary value of the anomalous 
quadrupole moment. In Fig. 1, I show q — (d — 6 — ah)/o? as a function of h for a* = 0.4 
(left panel) and a* = 1.4 (right panel). There are three distinct curves: one in the region 
b < a — M (curve A), one in the region a — M<b<a^M (curve B), and the last one 
for 6 > a + M (curve C). The equation g = {d — 6 — ab)/o? may have no real solutions for 
b (when |a*| < 1), and up to four distinct real solutions (for q sufficiently small). The Kerr 
metric is recovered for b^ — o? — . Two or more distinct solutions with the same M, a, 
and q correspond to objects with the same mass, spin, and quadrupole moment, but with 
different higher order moments. A particular solution requires prolate spheroidal coordinates 
if + 5 > and oblate spheroidal coordinates if + 5 < 0; otherwise, the constant k becomes 
an imaginary number. Generally speaking, slow-rotating solutions require prolate spheroidal 
coordinates, while fast-rotating solutions need oblate spheroidal coordinates; however, the 
critical value a^^^^ separating the two cases depends on q and can be either larger or smaller 
than 1. In the rest of the paper, I will study the evolution of the spin parameter of objects 
with q> —1 and therefore I will discuss the solutions of b belonging to the curve A and B. 
When I find two solution with the same quadrupole moment, I will call the solution with 
smaller MMSl, and the other, MMS2. 

3 Properties of the MMS solution 
3.1 Structure of the space-time 

The structure of the MMS space-time reminds one of the Tomimatsu-Sato and MN solu- 
tions [24, 48]. One finds naked singularities and closed time-like curves very close to the 
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massive object. The idea is that all these pathological features do not exist because inside 
the compact object, where the MMS metric, which is a vacuum solution, does not hold. 
Moreover, for |a*| < these pathological regions are always inside the inner radius of the 
disk and therefore they do not play any role in the discussion of the accretion process. 

The infinite redshift surface gu — determines the boundary of the ergoregion, which 
is not a pathological region and exists even around a rotating BH. Figs. 2, 3, and 4 show the 
ergosphere (red solid curve) for some MMS solutions with different value of the spin parameter 
and of the anomalous quadrupole moment. I use quasi-cylindrical coordinates pz, see App. A. 
The topology of the ergoregion can change significantly as a* and q vary. The same figures 
show also the boundary of the closed time-like curve regions (blue dotted curves), defined 
by ^(/)(/) = 0. For small deviations from the Kerr metric (Fig. 2, g = 0.01), the shape of the 
ergoregion is very similar to the Kerr case, see App. B, but there are one or two small holes, 
where gu < 0- Every hole of the ergoregion can be associated with a small region with closed 
time-like curves, even if the boundaries of the holes and of the closed time-like curve regions 
do not coincide. For g = 1.0 (Fig. 3), there are two disconnected ergoregions^ with no holes 
and one closed time-like curve region. For g = —1.0 (Fig. 4), one finds two ergoregions and one 
closed time-like curve region for the MMSl solutions (left panels) and one ergoregion with a 
hole and two distinct closed time-like curve regions for the MMS2 solutions (right panels). Let 
us notice that the MMSl solutions with (a*, q) = {(0.98, 0.01), (0.8, 1.0), (1.2, 1.0)} require 
prolate spheroidal coordinates, while all the other cases need oblate spheroidal coordinates. 
As the spin parameter increases, ergoregions and closed time-like curve regions move to larger 
p, but this is just an artifact of the coordinate system, see App. A. 

3.2 Geodesic motion 

The MMS space-time is stationary and axisymmetric. There are thus two constants of motion 
associated respectively with the t- and ^-coordinate; that is, the specific energy E and the 
specific axial component of the angular momentum L: 

E = -gut - gtcf)^ , L ^ gtcj^i ^ g^^cj) , (3.1) 

where the dot denotes the derivative with respect an affine parameter. From (3.1), we have 

j_ _ Egcfycf) + Lgtcf) ^ _ Egt^ + Lgu ^^^^ 
9t(f) ~ 9tt9(f)(i) ' 9tcf) ~ 9u9(t)ct) 

By substituting i and (j) into the equation of the conservation of the rest mass, g^yX^x^ = — 1, 
we find 

p27 / -2 -2 \ 27 

y (^ + r^j = y + = ^^ff ^' ^' ^) (3.3) 

where 

VeE = ^-^{L-ojEf-l. (3.4) 

^The fact that there are two disconnected ergoregions may depend on the coordinate system, as in prolate 
spheroidal coordinates the Schwarzschild radius r — M reduces to the segment \p\ < \/ — and z = 
and the region with r < M is not included, see App. A. The existence of these ergoregions should also be 
taken with caution, because close to the pathological region with closed time-like curves. 
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Since the left-hand side of Eq. (3.3) is non-negative, the motion of a test-particle is restricted 
to the region V^ff > 0. The zeros of the effective potential Vefr are shown in Figs. 5 and 6 
(orange dashed-dotted curves) for a* = 1.2, g = 1.0 (MMSl), and different values of E and 
L. They indicate the boundary of the allowed regions of the motion. The blue dotted curves 
denote the boundaries of the closed time-like curve regions and must be inside the object. 
In Fig. 5, top panels, oblate spheroidal coordinates xy are used. The bottom panels are 
enlargements of the region close to the compact object in quasi-cylindrical coordinates pz. 

For E = 1 and large L (right panels in Fig. 5), we find an unbound region extending 
to infinity and a plunging region close to the massive object. As L decreases, the unbound 
region moves closer to the plunging region and a new allowed region of motion appears 
(central panels of Fig. 5). The latter has no counterpart in the Kerr space-time, while it 
exists in the MN case [24]. For smaller L (left panels of Fig. 5), the unbound region merges 
with the plunging one. 

For E < 1, the situation is similar. In general, one find a bound region and a plunging 
region (central panel of Fig. 6). As L decreases, the bound region approaches the plunging 
region and eventually merges with the latter (left panel of Fig. 6). On the other hand, for 
larger values of L there is the plunging region only (right panel of Fig. 6). 

In conclusion, the behavior of the geodesic motion is not very different from the Kerr 
and MN space-times. In particular, one can expect that the standard picture of the accretion 
process is correct: the particles of the gas approach the compact objects by losing energy 
and angular momentum. When they are at the innermost stable circular orbit (ISCO), they 
quickly plunge to the object^. 

In addition to L, and the rest mass, in the Kerr space-time there is a fourth constant 
of motion, the Carter constant C [49]. Its existence allows for the full separability of the 
equations of motion and it is a major source of research in modern day EMRI modeling. For 
a generic space-time, the existence of an effective fourth constant of motion can be checked 
by studying the Poincare maps for geodesic motion of bound orbits: if they are all closed 
curves, there is a fourth constant of motion [24]. In analogy with the case of the MN solution 
discussed in [24], presumably even in the MMS space-time all orbits are strictly speaking 
chaotic, and no true fourth integral exists, but most of the orbits appear to be regular and 
one can find a quantity that is nearly invariant along them. 

4 Evolution of the spin parameter 

In this section, I compute the evolution of the spin parameter of the compact object due 
to the accretion process from a thin disk following Ref. [50]. I assume that the disk is 
on the equatorial plane of the object^ and that the disk's gas moves on nearly geodesic 
circular orbits. The gas falls to the central object by loosing energy and angular momentum. 

^This picture is surely correct for q > 0, as the radius of the ISCO is always determined by the orbital 
stability along the radial direction. For ^ < 0, the radius of the ISCO may depend on the stability of the 
orbit along the vertical direction, see e.g. the discussion in Refs. [33, 38]. In this second case, the gas particles 
reach the ISCO and then leave the equatorial plane. They plunge to the compact object after having lost 
additional energy and angular momentum. Such an additional loss of energy and angular momentum seems 
to be negligible in the description of the accretion process, but that should be checked in future studies by 
performing hydrodynamical and magnetohydrodynamical simulations. 

"^This assumption is correct for long term accretion onto a super-massive object at the center of a galaxy, 
since the alignment timescale of the spin of the object with the disk is typically much shorter than the accretion 
timescale [38]. 
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a* = 0.98,^ = 0.01 (MMSI) 



a. = 0.98,^ = 0.01 (MMS2) 




12 3 

a* = 1.2,^ = 0.01 (MMS1) 





12 3 

a* = 1.2,^ = 0.01 (MMS2) 




Figure 2. Space-time structure of the MMS solution for q = 0.01 and spin parameter a* = 0.98 
(top panels) and a* = 1.2 (bottom panels). The solid red curves denote the infinite redshift surface 
^tt = 0, defining the boundary of the ergoregion of the space-time. The dotted blue curves define the 
boundary of the causality violating region, where < 0. p and z are given in units of M = 1. 

a. = 0.8, ^=^.0 (MMSI ) a. = 0.8, g = 1 .0 (MMS2) 





a. = 1.2,g=1.0 (MMS2) 




Figure 3. As in Fig. 2, for ^ = 1.0 and spin parameter a* = 0.8 (top panels) and a* = 1.2 (bottom 
panels). 



When it reaches the ISCO, it plunges to the massive body. If the gas is "absorbed" by the 
compact object, with no further emission of radiation, the compact object changes its mass 
by 6M = EiscoSm and its spin angular momentum by 6 J = LiscoSm, where ^^isco ^ind 
^ISCO are respectively the specific energy and the specific angular momentum of the gas 
particle at the ISCO, while 6m is the gas rest-mass. The evolution of the spin parameter 
turns out to be governed by the following equation [50] 

dCi* 1 -^ISCO ry / . . N 

rflnM MEisco 

In the case of a body with a solid surface made of ordinary matter (i.e. protons, neutrons. 
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a* = 0.8,g = -1.0 (MMS1) 



a* = 0.8, g = -1.0 (MMS2) 





a* = 1.2, g = -1.0 (MMS1) 



12 3 

a* = 1.2.^1 = -1.0 (MMS2) 





E = 1.00, L = 3.1 



Figure 4. As in Fig. 3, for q = —1.0. 

E = 1.00, L = 3.6 



E = 1.00, L = 4.0 




Figure 5. Effective potential for geodesic motion around a super-spinning object with a* = 1.2 and 
q = 1.0 (solution MMSl) for E = 1.00 and L = 3.1 (left panels), L = 3.6 (central panels), and L = 4.0 
(right panels). Top panels: x?/-plane. Bottom panels: enlargement of the region closer to the object 
on the pz-plane. The dashed-dotted orange curves denote the zeros of the effective potential; dotted 
blue curves indicate the boundary of the causality violating region, where < 0. I/, p, and z are 
given in units of M = 1. 



and electrons), this picture is not rigorously correct: the gas particles hit the surface of the 
body and release their gravitational energy in form of radiation. They are also accumulated 
on the surface and may undergo nuclear reactions, with the production of bursts. In such a 
situation, the computation of the evolution of the spin parameter is much more complicated 
and model dependent. On the contrary, if the compact object is a BH, the gas particles are 
really absorbed without further emission of radiation, since the BH has no solid surface, and, 
once the particles are behind the event horizon, their radiation cannot escape to infinity. For 
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Figure 6. As in Fig. 5, for E = 0.94 and L = 2.8 (left panel), L = 2.9 (central panel), and L = 3.1 
(right panel) on the x?/-plane. 



a BH, one can analytically integrate Eq. (4.1) and obtain [50] 

forM < VqMq 



a* 



2 Mo 



M2 



4- a/18^ -2 



1 forM>\/6Mo, (4.2) 



assuming an initially non-rotating BH with mass Mq. The equilibrium is reached for — 1 
after the BH has increased its mass by a factor \/6 ^ 2.4. Including the effect of the radiation 
emitted by the disk and captured by the BH, one finds ^ 0.998 [34], because radiation with 
angular momentum opposite to the BH spin has larger capture cross section. The presence 
of magnetic fields in the plunging region may further reduce this value to ^ 0.95 [51, 52], 
by transporting angular momentum outward. 

For astrophysical BH candidates, we can presumably use Eq. (4.1). They indeed seems 
to be capable of absorbing the accreting gas without apparent loss of energy and angular 
momentum: no electromagnetic emission from their surface is observed, see e.g. Refs. [53- 
55]. For a generic stationary and axisymmetric space-time, Ei^qq and Lisco can be computed 
numerically, as described in [33] . Here I assume that g is a constant and does not change as 
the object increases its mass and is spun up, but this is not true in general. For example, 
in the case of neutron stars, q depends on the mass of the body [47]. If the astrophysical 
BH candidates were objects with an anomalous quadrupole moment depending strongly on 
the mass of the body, objects with different mass could have a different equilibrium spin 
parameter. 

Neglecting the effect of the radiation captured by the object and the presence of magnetic 
fields in the plunging region, in Fig. 7 I show da^/dlnM as a function of the spin parameter 
a* for some values of the anomalous quadrupole moment q. The value of the equilibrium 
spin parameter al^ for these cases is reported in Tab. 1. In Fig. 7, the solid red curve is for 
a BH (q = 0) and has been computed by using Boyer-Lindquist coordinates. Indeed, even if 
the Kerr metric is included in the MMS solution, neither the form of the metric presented 
in [44, 45], nor the one discussed in Sec. 2 are suitable for numerical calculations, as they 
would require that the parameter b is an imaginary number. For a similar reason, for a 
given q ^ we cannot study the evolution of the spin parameter from a* = to a* = a*^, 
but from some value a^^ > 0, depending on q. When da^/dlnM > 0, the accretion process 
spins the compact body up. When da^/dluM < 0, the compact body is spun down. The 
equilibrium spin parameter is thus the one for which da^/dln M = 0. In Fig. 8, I show the 
quantity da^/dlnM for a BH (green dotted curve) and for three compact objects with the 
same anomalous quadrupole moment q = ±1.0, but different higher order moments: the 
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q = -0.01 (MMS1) 
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Figure 7. da^/dlnM as a function of a* for different value of the anomalous quadrupole moment q. 
The red solid curve is for the case of a BH = 0). Left panel: cases with ^ > 0. Right panel: cases 
with q < 0. 



Kerr 
MN 
MMS1 
MMS2 





Figure 8. Comparison of the function da^/dliiM for the case of the MN solution with the MMSl 
and MMS2 solution with the same value of the anomalous quadrupole moment. Left panel: q = 1. 
Right panel: q = —1. The MN solution can be used till = 1. The green dotted curve is for the 
case of a BH. 



subclass of MN metrics discussed in Ref. [38] (red solid curve), MMSl (dark-blue dotted- 
dashed curve), and MMS2 (light-blue dotted-dashed curve). 

Fig. 9 shows the evolution of the spin parameter Qj^ as a function of M/Mq, the mass 
to initial mass ratio. Eq. (4.1) is numerically integrated from a^^, with M/Mq equal to the 
one of the Kerr solution with the same spin parameter. 

5 Discussion 

As the values of £^isco and Lisco in Eq. (4.1) depend on the metric of the space-time, 
compact objects with different quadrupole (or higher) moment have a different equilibrium 
spin parameter al^. Figs. 7, 9, and Tab. 1 can be qualitatively understood in term of the 
inner radius of the disk rin'. for a given spin parameter, da^/dluM depends on i^isco/^'iSCO? 
which typically increases/decreases if r^n is larger /smaller. Fig. 10 shows the value of the 
inner radius of the disk in Schwarzschild coordinates (see Eq. (A. 5) in App. A) as a function 
of the spin parameter. For instance, the MMSl solutions with ^ < have at^ < 1; as we 
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Figure 9. Evolution of the spin parameter a* as a function of M/Mq. The color and the style of 
the curves for any q is the same of Fig. 7 (the cases q = ±0.001 are not shown because for them the 
evolution of a* can be studied only in a very limited interval of the spin parameter). The red solid 
curve is for the case of a BH {q = 0). Left panel: cases with q > 0. Right panel: cases with q < 0. 



q 


af (MMSl) 


a7 (MMS2) 


af^ (Kerr) 


1.0 


1.44 


1.56 




0.1 


1.10 


1.17 




0.01 


1.02 


1.04 




0.0 






1.00 


-0.01 


0.97 


0.997 




-0.1 


0.92 


1.02 




-1.0 


0.85 


1.23 





Table 1. Equilibrium spin parameter al^ for the cases shown in Fig. 7. 

can see, their inner disk's radius is smaller than the one in Kerr space-time for |a*| < al^. 
The curves for the MMSl solutions q = —0.01 and —0.1 in the right panel of Fig. 10 stop at 
Tin = M, which is equivalent to x = 0, because our coordinates cannot describe the space- 
time at smaller radii (0 < x < +oc). These two cases have anyway to be taken with caution: 
there are regions with closed time-like curves with Schwarzschild radius larger than M, so 
even the surface of the compact object should be probably larger than M. Let us notice, 
however, that this happens for a* > al^ and it can be neglected in the study of the evolution 
of the spin. As shown in Fig. 10, for sufficiently high values of the spin parameter, the inner 
radius of the disk always increases as a* increases. This is not a surprise, because even in 
the Kerr space-time the radius of the ISCO reaches a minimum for a* ^ 1.089, and then it 
increases as a* increases [56]. 

The value of the inner radius of the disk depends inevitably on the choice of the co- 
ordinates; the comparison of this quantity for objects with different value of q is thus not 
so meaningful. From this point of view, a more interesting quantity is the angular velocity 
of a gas particle at the inner radius of the disk, which is shown in Fig. 11. As we can see, 
the angular frequency at r^n can be very high around a very fast-rotating BH, while it is 
significantly lower for another compact object, regardless of its spin parameter. This can be 
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an interesting observational feature to test the Kerr nature of astrophysical BH candidates, 
because it implies that it is very difficult for an object that is not a BH to mimic a very 
fast-rotating BH. In other words, among the stationary and axisymmetric space-times, the 
Kerr solution is a very special case with peculiar properties and this is particularly true for 
a* ^ 1. Generally speaking, the accretion process onto an object with non-Kerr quadrupole 
moment looks more like the accretion onto a slow-rotating BH, but it can hardly present 
some features of a fast-rotating BH. 

For fast-rotating objects, even deformations beyond the quadrupole moments are im- 
portant for the properties of the space-time. In Fig. 12, I show the inner radius of the disk 
and the corresponding angular frequency of a gas particle around objects with a* = 0.99 as 
a function of the anomalous quadrupole moment. The three curves represent three objects 
with same mass, spin, and quadrupole moment, but different higher order moments: the 
subclass of MN solutions discussed in [38], the MMSl solution, and the MMS2 solution. 

As already noticed in [38], Eq. (4.1) provides the correct evolution of the spin parameter 
if the compact object does not become unstable before reaching the equilibrium value . 
More in general, the accretion process can spin the body up, but there may exist other 
processes that spin it down. This is what should happen, for instance, in the case of a 
neutron star: the gas of the accretion disk can spin the neutron star up, but, when the 
rotational frequency of the object exceeds ^ 1 kHz, there are unstable modes that spin the 
neutron star down through the emission of gravitational waves [57]. If something similar 
happens for BH candidates, the maximum value of the spin parameter would be determined 
by the internal structure of these objects. The latter may spin down by emitting gravitational 
waves, presumably as a burst, potentially detectable by future experiments. 

The fact that the accretion process can spin a compact object up to a* > 1 can be 
relevant for the super-massive BH candidates at the center of galaxies, while it should be 
negligible for stellar-mass objects in X-ray binary systems. In general, the value of the spin 
parameter of a compact object is determined by the competition of three physical processes: 
the event creating the object, mergers, and gas accretion. For the stellar-mass BH candi- 
dates in X-ray binary systems, the value of their spin should reflect the one at the time of 
their creation. If they belong to low-mass X-ray binary systems, even swallowing the whole 
stellar companion they cannot change significantly their spin, because the mass of the stellar 
companion is much smaller. If they are in high-mass X-ray binary systems, even accreting at 
the Eddington limit they do not have enough time to grow before the explosion of the com- 
panion. On the contrary, for the super-massive objects in galactic nuclei the initial spin value 
is completely unimportant, as they have increased their mass by a few orders of magnitude 
from the original one. In the case of prolonged disk accretion, the object has the time to align 
itself with the disk and the process of gas accretion should dominate over mergers [36] . This 
picture is supported even by current estimates of the mean radiative efficiency of AGN [37] . 
As shown in [38], the alignment timescale is the same for BHs and objects with non-Kerr 
quadrupole moment, as long as \q\ <C 100. 

6 Conclusions 

The final product of the gravitational collapse of matter is thought to be a black hole and 
there are astrophysical evidences for the existence of dark objects that are too compact and 
too heavy to be relativistic stars or clusters of non-luminous bodies. In 4-dimensional general 
relativity, a black hole is completely specified by its mass M and by its spin parameter a*, 
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Figure 10. Inner radius of the accretion disk as a function of the spin parameter a* for different 
values of the anomalous quadrupole moment q. The color and the style of the curves for any q is the 
same of Fig. 7. Left panel: cases with g > 0. Right panel: cases with q < 0. Inner radius of the disk 
in Schwarzschild coordinates and in units of M = 1. 




Figure 11. Angular velocity at the inner radius of the accretion disk as a function of the spin 
parameter a* for different values of the anomalous quadrupole moment q. The color and the style 
of the curves for any q is the same of Fig. 7. Left panel: cases with q > 0. Right panel: cases with 
q < 0. Qin is given in units of M = 1. 

and it is subjected to the Kerr bound |a>,<| < 1, which is the condition for the existence of 
the event horizon. The accretion process can spin a black hole up to a* ^ 0.998 and at least 
some of the super- massive objects at the center of galaxies may be fast-rotating black holes 
with spin parameter close to this value. 

The Kerr black hole paradigm is still based on a set of unproven assumptions. They 
sounded reasonable forty years ago, but they are more questionable today. There is no di- 
rect evidence that black hole candidates have an event horizon, while there are theoretical 
arguments suggesting new physics appearing at macroscopic scales [7-12]. If the astrophys- 
ical black hole candidates are not the objects predicted by general relativity, they are not 
subjected to the Kerr bound. Interestingly, the accretion process onto a body with non-Kerr 
quadrupole moment typically spins the object up to a* > 1. 

In Ref. [38], I used the Manko-Novikov solution to describe the exterior gravitational 
field of a generic compact body. I showed that in most cases the equilibrium spin parameter is 
larger than 1. However, the Manko-Novikov solution is valid only for la^l < 1 and therefore 





Figure 12. Left panel: inner radius of the disk for bodies with spin parameter a* = 0.99 as a 
function of the anomalous quadrupole moment q for the MN (red solid curve), MMSl (blue dashed- 
dotted curve), and MMS2 (green dotted curve) solutions. Right panel: as in the left panel, for the 
angular velocity at the inner radius of the disk. Inner radius of the disk in Schwarzschild coordinates 
and in units of M = 1. is given in units of M = 1. 



it was impossible to study the accretion process for a* > 1 and figure out the properties 
of the space-time around a super-spinning body. In the present paper, I considered the 
Manko-Mielke-Sanabria- Gomez solution, which can describe the gravitational field outside 
a compact body with non-Kerr quadrupole moment and spin parameter either smaller and 
larger than 1. I studied the basic properties of the space-times when a* > 1, I discussed the 
evolution of the spin parameter, and I found its equilibrium value, see Figs. 7, 9, and Tab. 1. 
For fast-rotating objects, the accretion processes onto a black hole and onto a generic body 
are quite different. For example, only around a black hole can the inner radius of the disk be 
very small, while the angular frequency of the gas particles at the innermost circular orbit 
can be significantly higher than the one around another object. 

The fact that the accretion process in a non-Kerr background can spin the compact body 
up to a* > 1 is relevant for the research devoted to figure out how future observations can 
test the Kerr nature of the current astrophysical black hole candidates. The possibility that 
these objects can have spin parameter larger than 1 cannot be ignored and this is particularly 
true for experiments like LISA, whose detection relies on matched filtering. 
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A Coordinate systems 

The canonical form of the line element of a generic stationary and axisymmetric space-time 
in quasi-cylindrical coordinates pz is 

27 2 
ds"^ = -/ {dt - ud(j)f + ^ (V + dz^) + _ (^A.l) 
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The relation between quasi-cylindrical coordinates pz and prolate spheroidal coordinates xy 
is qiven by 



p = ^2)^ z^kxy, (A.2) 

where is a constant. The inverse relation is 

X = — ^- , y = , A.3) 

where R± — ^y p"^ + {z ± k^. Through the formal transformation X y ix and k —ik, 
where i is the imaginary unit, one changes the prolate spheroidal coordinates into oblate 
spheroidal coordinates. The relation between the quasi-cylindrical coordinates and the oblate 
spheroidal coordinates is therefore given by 

p k\/{x^ + 1)(1 - y^) , z^kxy. (A. 4) 

Lastly, the relation between quasi-cylindrical coordinates and Schwarzschild coordinates 
is [58] 



p = Vr2^^2MrTa2sin6>, z = (M - r) cos 6> , (A. 5) 

where M is the mass and a — J/M the specific spin angular momentum. 

Prolate spheroidal coordinates are suitable for describing the space-time around slow- 
rotating objects, while oblate spheroidal coordinates are suitable in the case of fast-rotating 
objects. In the special case of the Kerr space-time, prolate spheroidal coordinates can be used 
only for BHs, oblate spheroidal coordinates only for Kerr naked singularity (see App. B). The 
left (right) panel of Fig. 13 shows some curves with constant prolate (oblate) coordinate x on 
the pz-plane for k = 1.5 M. Fig. 14 shows instead some curves with constant Schwarzschild 
radial coordinate, still on the pz-plane, for a — 0.8 M (left panel) and a — 1.2 M (right panel). 
In the former case, the quasi-cylindrical coordinates cover the space with Schwarzschild 
radial coordinate r > r/f = M + \/ — o? and the surface with radial coordinate r — th 
corresponds to the segment p = and \z\ < \JM^ — a^. Here tu is equal to the Schwarzschild 
radial coordinate of the even horizon of a Kerr BH with spin a. When \a\ > M, the quasi- 
cylindrical coordinates can be used to describe the region r > M (otherwise there is not a 
one-to-one correspondence between the two coordinate systems) and the circle with radius 
r — M reduces to the segment \p\ < \J cv^ — and z = on the pz-plane. For non-Kerr 
objects, prolate (oblate) spheroidal coordinates are still adequate only for slow-rotating (fast- 
rotating) objects, but in general the value of the spin a separating the two cases is not M 
any more. 

B Kerr space-time in spheroidal coordinates 

In prolate spheroidal coordinates, the line element (A.l) is 



,,,0 fc^e^T .9 9^ / dx^ dir \ 

-f (dt - cod<pf + — - y') (^231 + j + 

+ ^{x'-l){l-y')d<p\ (B.l) 
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Figure 13. Curves with constant prolate spheroidal coordinate x (left panel) and oblate spheroidal 
coordinate x (right panel) on the p2:-plane. p and z are given in units of M = 1. k = 1.5 M. 



= 0.8 



: 1.2 



N 




N 




-0.5 



Figure 14. Curves with constant Schwarzschild radial coordinate on the /)2;-plane, for — 0.8 (left 
panel) and a* = 1.2 (right panel), p and z are given in units of M = 1. 



In the Kerr space-time 

where A, and C are given by 

A = k^{x^ - 1) - a^{l -y^), B = {kx + Mf + a^y^ , C = 2aM{kx + M) . (B.3) 

Here M is the mass, a — J /M \s the specific spin angular momentum, and k — \J — o?. 
Prolate spheroidal coordinates can be used to describe only BHs (|a>,<| < 1) and the surface 
X = 1 is the BH event horizon. The transformation x ^ ix and k —ik changes the prolate 
spheroidal coordinates into oblate spheroidal coordinates, which can be used to describe the 
space region around a Kerr naked singularity (|a>,<| > 1) with Schwarzschild radial coordinate 
r > M. Fig. 15 shows the infinite redshift surface gu = (red solid curve) and the event 
horizon (blue dotted curve) of a BH with a* = 0.8 in prolate spheroidal coordinates xy 
(left panel), quasi-cylindrical coordinates pz (central pane), and Schwarzschild coordinates 
r0 (right panel). Fig. 16 shows the infinite redshift surface of a Kerr naked singularity with 
a* = 1.2 in oblate spheroidal coordinates xy (left panel), quasi-cylindrical coordinates pz 
(central pane), and Schwarzschild coordinates r0 (right panel). 
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Figure 15. Infinite redstiift surface gu = (red solid curves) and event tiorizon (blue dotted curves) 
in Kerr space-time with a* = 0.8 in prolate spheroidal coordinates xy (left panel), quasi-cylindrical 
coordinates pz (central panel), and Schwarzschild coordinates r6 (right panel), p, and r are given 
in units of M = 1. 
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Figure 16. Infinite redshift surface gu = in Kerr space-time with a* = 1.2 in oblate spheroidal coor- 
dinates xy (left panel), quasi-cylindrical coordinates pz (central panel), and Schwarzschild coordinates 
rO (right panel), p, z, and r are given in units of M = 1. 
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